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Abstract  —  A  simulation  and  experimental  study  has  been 
conducted  on  equivalent  dipole  layer  imaging  (EDLI)  of  brain 
electric  sources  from  EEG.  Using  the  three-sphere 
inhomogeneous  head  model,  the  performance  of  the  EDLI  was 
rigorously  evaluated  for  a  variety  of  brain  source 

configurations  under  different  noise  levels.  The  present 
simulation  results  demonstrate  the  excellent  performance  of 
the  EDLI  in  mapping  and  imaging  the  underlying  cortical 
sources  with  much  enhanced  spatial  resolution,  as  compared  to 
the  scalp  potentials.  Human  experiments  were  further 

conducted  to  examine  the  feasibility  of  EDLI.  Pattern  reversal 
visual  evoked  potentials  (VEP)  were  recorded  from  94 
electrodes  and  the  brain  electric  sources  at  P100  were 
estimated.  The  VEP  experiments  demonstrate  that  the  present 
EDLI  can  eliminate  the  misleading  far  field  in  the  scalp 
potential  map,  localize  and  map  the  underlying  cortical  sources 
induced  by  the  visual  stimuli. 
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I.  Introduction 

Although  conventional  EEG  offers  excellent  temporal 
resolution  in  resolving  rapidly  changing  pattern  of  brain 
electric  activities,  its  spatial  resolution  is  limited  by  the 
smoothing  effect  of  the  head  volume  conductor,  especially 
due  to  the  very  low  conductivity  skull  layer  [1], 

A  number  of  investigators  have  attempted  to  improve  the 
spatial  resolution  of  the  EEG  by  solving  the  EEG  inverse 
problem  [2-13].  Of  particular  interest  is  the  recently 
developed  cortical  imaging  technique  (CIT)  [7-13],  in  which 
an  explicit  biophysical  model  of  the  passive  conducting 
properties  of  a  head  is  used  to  deconvolve  a  measured  scalp 
potential  distribution  into  a  potential  distribution  on  the 
cortical  surface.  The  CIT  needs  to  construct  a  virtual  dipole 
layer  (DL)  inside  the  cortical  surface  and  a  closed  DL 
equivalently  represents  all  the  brain  electric  sources 
bounded  by  its  surface  [10-13],  The  inverse  procedure 
estimates  the  equivalent  DL  from  the  scalp  EEG,  then  the 
cortical  potentials  are  reconstructed  by  solving  the  forward 
problem,  from  the  estimated  equivalent  DL  to  the  cortical 
potentials  [7-13]. 

Although  as  an  intermediate  result  of  CIT,  the  estimate 
of  the  DL  may  also  facilitate  the  equivalent  representation 
of  brain  electric  sources,  previous  CIT  approaches  had  to 
further  calculate  the  cortical  potentials,  due  to  the  lack  of  the 
forward  theory  on  the  equivalent  DL.  Recently,  we  have 
proposed  the  equivalent  dipole  layer  imaging  (EDLI),  which 
provides,  for  the  first  time,  the  forward  theory  of  the 
equivalent  DL.  Based  on  this  theory,  the  brain  electric 


sources  can  be  directly  imaged  on  the  equivalent  DL  from 
EEG  measurements,  thus  eliminating  the  need  for  the 
forward  calculation  of  cortical  potentials  in  the  CIT.  In  the 
present  study,  we  conducted  both  computer  simulations  and 
human  experiments  to  rigorously  evaluate  the  performance 
of  EDLI  and  its  feasibility  to  image  brain  activation. 

II.  METHODS 

A.  Forward  solution  of  the  equivalent  DL 

Considering  an  arbitrary-shaped  head  model  as  illustrated 
in  Fig.  1,  where  S’  represents  the  scalp  surface  and  S 
represents  the  virtual  DL.  Denote  V0  as  the  volume  bounded 
within  S,  and  V  the  volume  bounded  by  S  and  S’.  If  all 
brain  electric  sources  are  in  V0  and  there  is  no  active  source 
in  V,  then  for  any  given  point  x  within  V,  the  potential 
function  (p(x)  can  be  expressed  as  [14]: 

tp(x)  =  f  [G{x,x'){o  -)-(<n (p(x'))^G^,X  ^ds'  (1) 

Js  on'  on' 

where  x’  is  on  S,  d/dri  is  the  normal  derivative  at  surface  S, 
a  is  the  conductivity  in  V0,  G(x,x’)  is  the  Green  function 
which  is  the  potential  solution  of  a  point  current  source 
inside  a  volume  conductor  model. 


Fig.  1  Illustration  of  DL  in  an  arbitrary-shaped  head  model 

Different  from  the  conventional  Green  function  method 
[14],  Eq.  (1)  can  be  solved  by  keeping  G  unchanged,  while 
modifying  function  cp  by  decomposing  it  into  two  potential 
components  produced  by  the  sources  in  volumes  V  and  V0, 
respectively.  Since  the  sources  in  volume  V  make  no 
contribution  to  the  surface  integral  in  Eq.  (1),  we  can 
replace  cp(x)  with  a  proper  (p’(x)  so  that: 

(p'(x')  =  (p'(x')v  +(p(x')vo,  for*’  on  S  (2a) 

dcp'(x')  _  dcp'(x')v  |  dcp(x')v „  _Q  forx’  on  iS  (2b) 

dn'  dn'  dn' 

Then  Eq.  (1)  can  be  simplified  to: 

cp(x)  =  |  (<7<p'(x'))  _  f  (atp' (x'))Gd(x, x')ds' * 

Js  dn'  Js 

for  a  ’  oil  S  (3) 

where  Gd  =-dG/dn'  is  the  Green  function  of  a  dipole  on 
surface  S  oriented  normally  outward.  Therefore,  Eq.  (3)  can 
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be  considered  as  an  equivalent  DL  model,  with  equivalent 
dipole  source  density  of  f  =  o(p'  (x  ) ,  and  with  Green 

function  of  G,t.  Eq.  (3)  is  valid  for  an  arbitrary  geometric 
model  and  numerical  methods  can  be  applied  to  calculate 
G,t.  Specially,  when  the  equivalent  DL  is  considered  to  be  a 
spherical  surface,  the  closed  solution  of  f.  can  be  obtained, 
with  (p'(x')  be  the  surface  potential  produced  by  a  dipole 
inside  a  homogeneous  conducting  sphere  with  boundary 
condition  as  described  in  Eq.  (2a,b).  This  implies  that  fd  is 

proportional  to  the  potential  over  the  same  spherical  surface 
when  the  exterior  space  of  the  DL  is  replaced  by  air.  In 
another  word,  the  equivalent  DL  source  density  distribution 
may  reflect  the  potential  distribution  over  this  layer,  had  the 
upper  medium  being  removed  during  the  open-skull 
surgery. 

B.  Inverse  problem  ofEDLI 

The  proposed  EDLI  constructs  the  equivalent  dipole 
source  distribution  over  the  DL,  which  is  essentially  not 
affected  by  the  low-conductivity  skull  layer,  and  has  higher 
spatial  resolution  as  compared  to  the  scalp  potential  map. 
Discretize  Eq.  (3)  as: 

<p(x)  =  £  (/,  (Ms,  )Gd  (. x ,  x ' )  ■  for  *’  on  S  (4) 

Denote  the  source  density  weighed  by  discrete  grid  area 
f,  ( i)dsj  as  the  equivalent  DL  source  strength.  Then  Eq.  (4) 

linearly  relates  the  potential  (p  for  x  in  V,  with  the  equivalent 
DL  source  strength,  by  the  discrete  Green  function  Gd. 
Specially,  when  (p  is  the  measured  scalp  potential,  Eq.  (8) 
can  be  written  in  matrix  form: 

=  AF  (5) 

where  (j)s  is  the  vector  of  measured  scalp  potentials,  F  is 

the  vector  of  equivalent  DL  source  strength,  A  is  the  lead 
field  matrix  obtained  by  evaluating  Gd. 

Denote  the  pseudo-inverse  of  A  as  A+ ,  the  equivalent  DL 
source  strength  can  be  estimated  as: 

F  =  A+$s  (6) 

The  proposed  EDLI  is  based  on  the  inverse  estimation  of 
vector  F ,  and  the  truncated  singular  value  decomposition 
(TSVD)  was  applied  to  improve  the  numerical  stability  of 
the  inverse  problem  [15]. 

C.  Simulation  and  experimental  protocols 

Computer  simulations  based  on  the  three-sphere 

inhomogeneous  head  model  [16]  were  conducted  to  evaluate 
the  EDLI.  The  normalized  radii  of  the  equivalent  DL,  the 
brain,  the  skull  and  the  scalp  spheres  were  taken  as  0.80, 
0.87,  0.92  and  1.0,  respectively.  The  normalized 

conductivity  of  the  scalp  and  the  brain  was  taken  as  1.0,  and 
that  of  the  skull  as  0.0125.  Lour  unit-moment  dipole  sources 
(radial  or  tangential)  with  varying  eccentricity  were  used  to 
represent  four  well-localized  brain  electric  activity.  The 
forward  solution  of  the  equivalent  DL  was  validated  by 
evaluating  the  relative  error  (RE)  and  correlation  coefficient 


(CC)  between  the  analytic  cortical  potentials  generated  by 
the  dipoles  and  the  cortical  potentials  produced  by  the 
equivalent  DL.  To  evaluate  the  inverse  solution  of  EDLI, 
Gaussion  white  noise  (GWN)  was  added  to  the  scalp 
potentials  generated  by  the  simulating  dipoles  to  simulate 
noise-contaminated  scalp  potential  measurement.  The 
source  strength  of  the  equivalent  DL  was  estimated  using 
Eq.  (6)  and  compared  with  the  forward  solution.  The  TSVD 
was  applied  to  solve  the  inverse  problem,  and  the  minimal 
product  (MINP)  method  was  used  to  determine  the 
truncation  parameter  in  TSVD  [12]. 

Human  visual  evoked  potentials  (VEPs)  were  recorded 
from  two  healthy  subjects  according  to  a  protocol  approved 
by  UIC/IRB.  Visual  stimuli  were  generated  by  the  STIM 
system  (Neuro  Scan  Labs).  94-channel  VEP  signals 
referenced  to  right  earlobe  were  amplified,  band  pass 
filtered  and  acquired  by  the  ESI  systems  (Neuro  Scan  Labs). 
The  electrode  locations  were  measured  using  Polhemus 
Lastrack  (Polhemus  Inc.)  and  best  fitted  on  the  spherical 
surface  with  unit  radius.  Half  visual  field  pattern  reversal 
checkerboards  with  reversal  interval  of  0.5  sec  served  as 
visual  stimuli  and  400  reversals  were  grand  averaged  to  get 
VEP  signals.  The  EDLI  maps  at  PI 00  were  estimated  and 
compared  to  the  corresponding  scalp  potential  maps. 

III.  Results 

A.  Evaluation  of  the  forward  solution  of  equivalent  DL 

Tab.l  shows  the  effect  of  source  eccentricity  on  the 
accuracy  of  forward  solution  of  equivalent  DL,  which  was 
discretized  into  1280  grids.  Pour  radial  or  tangential  dipoles 
atr-(±sin(;r/6),0,cos(;r/6))  and  r •  (0,±sin(7r/6), cos(;r/6)) 
were  used  to  simulate  brain  electric  sources,  where  r  is  the 
eccentricity  of  the  dipoles  and  ranges  from  0.30  to  0.75. 
Tab.  1  clearly  indicated  that  for  both  radial  and  tangential 
dipoles,  very  low  RE  and  high  CC  values  were  achieved  for 
most  source  eccentricities,  except  for  r  =  0.75  when  sources 
are  very  close  to  the  DL. 


Tab.  1  Effect  of  source  eccentricity  (ECC) 


ECC 

0.20 

0.30 

0.40 

0.50 

0.60 

0.70 

0.75 

RE-rad 

0.0281 

0.0281 

0.0282 

0.0290 

0.0336 

0.0652 

0.1885 

CC-rad 

0.9996 

0.9996 

0.9996 

0.9996 

0.9994 

0.9979 

0.9842 

RE-tag 

0.0279 

0.0278 

0.0280 

0.0285 

0.0300 

0.0463 

0.1759 

CC-tag 

0.9996 

0.9996 

0.9996 

0.9996 

0.9996 

0.9989 

0.9844 

Tab.  2  shows  the  effect  of  the  grid  density  on  the 
accuracy  of  forward  solution  of  equivalent  DL.  Pour  radial 
or  tangential  dipoles  (as  described  above)  were  fixed  at 
t-0.60,  and  the  equivalent  DL  was  constructed  with  varying 
grid  densities  (320,  640,  1280,  2560,  and  5120).  Tab.  2 
clearly  revealed  that  for  both  radial  and  tangential  dipole 
sources,  the  denser  grid  of  the  equivalent  DL,  the  lower  RE 
and  higher  CC  values  can  be  achieved. 

B.  Imaging  brain  electric  sources  using  EDLI 

Pig.  2  shows  a  typical  example  of  inverse  estimation  of 
simulated  brain  electric  sources  using  the  EDLI  approach. 
Pig.  2(a)  displays  the  noise  (5%  GWN)  contaminated  scalp 
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Tab.  2  Effect  of  grid  density  of  the  equivalent  DL 


Grid 

320 

640 

1280 

2560 

5120 

RE-rad 

0.2685 

0.2009 

0.0336 

0.0294 

0.0020 

CC-rad 

0.9646 

0.9800 

0.9994 

0.9996 

1.0000 

RE-tag 

0.2596 

0.1857 

0.0300 

0.0243 

0.0019 

CC-tag 

0.9670 

0.9829 

0.9996 

0.9997 

1.0000 

potential 

map 

for  4 

radial 

dipoles 

located 

0.7  •  (±sin(;r /7),0, cos(;r /7))  and  0.7-(±sin(jr/7),0,cos(;r/7)). 
Notably,  the  scalp  potential  map  was  severely  blurred  and 
distorted  by  the  head  volume  conductor  and  additive  noise. 
The  forward  solution  of  the  EDLI  and  the  inversely 
estimated  EDLI  maps  were  respectively  shown  in  Fig.  2(b) 
and  (c).  Without  the  blurring  effect  caused  by  the  skull 
layer,  the  EDLI  maps  clearly  revealed  the  underlying  dipole 
sources  with  much  enhanced  spatial  resolution  as  compared 
to  the  scalp  potential  map,  and  the  estimated  EDLI  map 
corresponds  well  to  the  forward  solution  of  EDLI  map. 


(a)  (b)  (c) 

Fig.  2  An  example  of  estimation  of  simulated  brain  electric  sources 
using  EDLI.  All  maps  are  normalized  for  display. 


Fig.  3  shows  the  application  of  EDLI  in  VEP  data 
analysis.  Subject  A  and  B  were  respecitively  given  left  and 
right  visual  field  stimuli.  Their  normalized  scalp  potential 
maps  (top-back  view)  at  PI 00  were  respectively  shown  in 
Fig.  3(a)  and  (b),  and  the  estimated  EDLI  maps  (top-back 
view)  were  respectively  shown  in  Fig.  3(c)  and  (d).  Notably, 
in  response  to  the  left  visual  field  stimuli,  a  dominant 
positive  potential  component  was  elicited  with  a  widespread 
distribution  on  the  left  scalp  (Fig.  3(a)).  While  in  response 
to  the  right  visual  field  stimuli,  a  widely  distributed  positive 
potential  component  was  dominant  on  the  right  scalp  (Fig.  3 
(b)).  On  the  other  hand,  the  estimated  EDLI  map  in  response 
to  the  left  visual  field  stimuli  revealed  a  dominant  and  more 
localized  activity  in  the  right  visual  cortex  (Fig.  3(c)). 
Consistently,  the  estimated  EDLI  map  in  response  to  the 
right  visual  field  stimuli  indicated  a  dominant  and  more 
localized  activity  in  the  left  visual  cortex  (Fig.  3  (d)). 


(c)  (d) 

Fig.  3  Application  of  the  EDLI  to  VEP  data  analysis.  All  maps  are 
normalized  for  display. 

IV.  DISCUSSION 

From  the  point  of  view  of  electromagnetic  theory,  a 
closed  surface  DL  can  equivalently  represent  all  the  electric 
sources  enclosed  by  its  surface.  By  constmcting  an 
equivalent  DL,  the  potential  everywhere  between  this  DL 
and  the  scalp  can  be  determined,  including  the  cortical 
potentials  as  in  CIT  [7-13].  We  have  recently  developed  the 
forward  theory  of  the  EDLI,  and  demonstrated  that  the 
equivalent  DL  source  density  is  proportional  to  the  potential 
over  the  same  surface  had  the  upper  medium  being 
removed.  This  has  made  it  possible  to  evaluate  the 
equivalent  source  distribution  directly  over  the  equivalent 
DL,  without  a  further  calculation  of  the  cortical  potentials  as 
in  CIT. 

In  the  present  study,  using  a  3-sphere  inhomogeneous 
head  model,  we  have  evaluated  the  forward  solution  of  the 
EDLI  by  computer  simulations.  The  results  indicate  that  the 
equivalent  DL  can  well-represent  the  simulating  dipole 
sources  with  different  eccentricity  and  with  radial/tangential 
orientation  (Tab.  1),  and  denser  grid  of  the  equivalent  DL 
yields  more  accurate  representation  of  the  brain  electric 
sources  (Tab.  2).  We  also  evaluated  the  EDLI-based  inverse 
problem  by  both  computer  simulations  and  experimental 
studies.  As  shown  in  Fig.  2,  the  estimated  EDLI  map 
corresponded  well  to  the  analytic  EDLI  map,  and  both  of 
them  clearly  resolved  the  underlying  dipole  sources  with 
much  enhanced  spatial  resolution  as  compared  to  the  noise 
contaminated  scalp  potential  map.  Promising  results  were 
also  obtained  from  EDLI  analysis  of  P100  of  VEP  data  (Fig. 
3).  It  is  widely  accepted  that  the  half  visual  field  stimuli 
activate  the  visual  cortex  on  the  contralateral  hemisphere  of 
the  brain.  But  paradoxically,  the  half  visual  field  stimuli 
elicited  stronger  positive  potential  distribution  over  the 
ipsilateral  side  of  the  scalp  (Fig.  3(a)-(b)),  which  might  be 
misinterpreted  as  ipsilateral  visual  cortex  activation  [17]. 
However,  the  estimated  EDLI  maps  clearly  indicated  that 
the  contralateral  visual  cortex  was  activated  (Fig.  3(c)-(d)), 
thus  effectively  eliminated  the  misleading  far  field  observed 
in  the  scalp  potential. 

In  summary,  the  promising  results  of  the  present 
investigation  suggest  the  EDLI  approach  may  offer  a 
promising  alternative  means  of  mapping  spatially 
distributed  brain  electric  activity  noninvasively.  Further 
improvement  of  the  EDLI  can  be  achieved  by  taking  into 
account  the  realistic  geometry  head  model  and  arbitrarily 
shaped  equivalent  DL. 
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